Resolving magnetic frustration in a pyrochlore lattice 
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CeFe2 is a geometrically frustrated ferromagnet that lies close to an instability at which a subtle 
change in the lattice symmetry couples to a transition to antiferromagnetism. We use x-ray diffrac- 
tion, diamond-anvil-cell techniques, and numerical simulation to identify the ground states and to 
quantitatively illustrate effects of competing magnetic energy scales and geometrical frustration on 
the magnetic phase diagram. Comparison of phase transitions under both chemical substitution and 
applied pressure suggests a general solution to the physics of pyrochlore rare earth inter-metallic 
magnets. 



PACS numbers: 75.30.Kz, 75.50.Ee, 75.20.En, 61.05.cp 

In a geometrically frustrated magnet not all pairs of 
magnetic ions can simultaneously assume their lowest 
energy configuration, and the system often has to choose 
among a large number of competing configurations. From 
a conceptual standpoint, frustrated systems are inter- 
esting because ostensibly minor phenomena can play an 
outsized role in dictating the overall magnetic structure. 
CeFe2 is a frustrated rare earth inter-metallic magnet 
with Laves crystal symmetry. Laves phase materials be- 
long to a broader class of pyrochlore lattice systems, 
which attract much interest as canonical examples of geo- 
metrically frustrated magnetism in three dimensions pQ. 
Many Laves magnets of both ferro- and antiferromag- 
netic symmetry are known [2]-[7] and, as we demonstrate 
here, the ease with which the magnetic symmetry can be 
switched is a direct result of geometrical frustration. In 
this work we quantitatively illustrate the effects of geo- 
metrical frustration on the generic magnetic phase dia- 
gram of CeFe2. 

The cubic Laves structure features transition metal 
ions on a pyrochlore lattice of corner-sharing tetrahcdra 
woven through by corrugated planes of rare earth ions, 
as illustrated in Fig. [I] Stoichiometric CeFe2 [8HTU] is 
ferromagnetic but undergoes a phase transition to an an- 
tifcrromagnet under light chemical doping by for example 
Al, Co, Ru or Ir [TTI - flrJj . Previous studies on stoichio- 
metric CeFe2 have found evidence for a high pressure 
phase, which by comparison with the doped systems was 
posited to be antiferromagnetic . Further support for 
this identification comes from the observation of antifer- 
romagnetic fluctuations within the ferromagnetic phase 
[SI [T7] . We show in Fig. [I] that the existing data for 
chemical doping and applied pressure can be collapsed 
onto a generic phase diagram, wherein both the ferromag- 
netic and the antiferromagnetic transition temperatures 
are plotted as a function of an effective tuning parame- 
ter rj. Interestingly, r\ does not vary monotonically with 



the lattice constant, and the antiferromagnetic phase can 
be found for materials with lattice constants both larger 
and smaller than the parent compound CeFe2. This ap- 
parent insensitivity to the lattice constant is at variance 
with the conventional notion that chemical substitution 
influences magnetism by exerting a positive or a negative 
"chemical pressure" , and suggests that the phase diagram 
is instead controlled by symmetry and the geometry of 
magnetic frustration. 

With sensitivity to both lattice and magnetic symme- 
try, and compatibility with high pressure diamond anvil 
cells [IH1 [H], synchrotron x-ray diffraction is the ideal 
technique to study the phase diagram of a geometrically 
frustrated antifcrromagnet (Supplemental Material). We 
plot in Fig.[2]the response of the crystal lattice to applied 
pressure P at low temperature. Near P = 1.5 GPa the 
lattice undergoes a transition in which it is compressed 
along one of the cubic body diagonals, resulting in a split- 
ting of the cubic crystal into four types of rhombohedral 
domains. Using high resolution x-ray diffraction it is pos- 
sible to index diffraction peaks by their rhombohedral 
domain type. The compressibility Bq is determined by 
fitting the lattice constant a(P) to a one parameter Birch 
equation [19], with B = 90 ± 4 GPa and 105 ± 5 GPa 
in the low and high pressure phases respectively. In the 
high pressure phase the unit cell angle a deviates from 
90° by 0.327 ± 0.002°. The data also show clear evi- 
dence for a regime of phase coexistence between 1.3 and 
1.8 GPa. The presence of a structural phase transition 
and a regime of phase coexistence are consistent with 
magnetic susceptibility studies of CeFe2 under pressure, 
which found the ferromagnetic phase boundary to be dif- 
ficult to pin down at low temperature [IT] . Note also that 
a similar rhombohedral distortion is known to accompany 
the ferromagnetic-to-antiferromagnetic phase transition 
in Ru-, A1-, and Co-doped CeFe2, with a in the range 
90.2° - 90.31° 0OUJ, further supporting the notion that 
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FIG. 1: (color online). Crystal structure and magnetic phase 
diagram of CeFe2. (a) Layered view of CeFe2 along (111). 
The structure consists of stacked sheets of 3e-Fe sites in a 
kagome lattice, alternating with corrugated sheets of Ce sites 
and lb-Fe sites in a triangular lattice. For clarity, the (111) 
axis is elongated and not all atoms are shown, (b-d) The 
twelve nearest neighbors for each site, (e) Generic magnetic 
phase diagram for CeFe2 showing the paramagnetic (PM), 
ferromagnetic (FM) , and antiferromagnetic ( AF) phases. The 
tuning parameter r\ is a linear combination of pressure P (in 
GPa) and doping x (in %) according to r\ — P + Ax, where 
A is a dopant-specific scaling factor: A — 20, 65, 55, and 55 
GPa/% for Co, Al, Ru, and Ir, respectively. All four dopants 
would individually preserve the Laves phase when fully re- 
place Fe in CeFe2 [9]. For pure CeFe2 under pressure and 
Co-doped CeFe2, the lattice constant shrinks with increasing 
77 |HH13ll26] , In contrast, the lattice expands for doping with 
Al, Ru and Ir pHTfJj . 
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FIG. 2: (color online). Evolution of the CeFe2 crystal lattice 
with pressure at 3.5 K. (a) Longitudinal (0 — 29) scans of 
the (220) and (311) Bragg peaks at three pressures: in the 
low-pressure cubic phase, in the coexistence regime, and in 
the high pressure rhombohedral phase. The peak splitting 
at high pressure is evidence of the rhombohedral distortion. 
At 1.8 GPa, the splitting of the (220) peak indicates phase 
coexistence, while the cubic (311) peak is indistinguishable 
from the rhombohedral (311) peak within our measurement 
resolution, (b) Dependence of the lattice constant a and the 
cell-axis angle a on pressure. The shaded area marks the 
phase coexistence regime. 



the magnetic phase diagram is linked to symmetry con- 
siderations rather than controlled by effective chemical 
pressure or the effects of disorder. 

Using non-resonant magnetic x-ray diffraction we 
searched for and found the high pressure antiferromag- 
netic phase in the form of (1/2, 1/2, l/2)-type Bragg 
diffraction peaks associated with antiferromagnetic pe- 
riod doubling 20 . The longitudinal magnetic peak 
widths are limited by the instrument resolution (Fig. [3^,) , 
meaning that the coherence length of the antiferromag- 
netic domains is at least 1500 A. This direct observation 
of antiferromagnetic order in compressed, stoichiometric 
CeFe2, not just in its doped analogues, is strong justifi- 
cation for the generic phase diagram drawn in Fig. [T] 

We show in Fig. [3j3 antiferromagnetic reflections at 
nine positions in reciprocal space, all corresponding to a 
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FIG. 3: (color online). Direct observation of antiferromag- 
netism at 3.3 GPa and 3.5 K. (a) Longitudinal scans of mag- 
netic peaks (red) normalized to the lattice peaks (blue) re- 
spectively, (b) Quantity |S/ m e lqr /2 ■ §2\ 2 calculated using 
Eq. for six measured diffraction peaks (black circles). Also 
shown are the sensitivity limits of null measurements at three 
positions (black error bars with downward arrows). Calcu- 
lated values (purple line) are given for the model described in 
the text. 



single rhombohedral domain at P = 3.3 GPa and T = 3.5 
K. The nine positions yielded six measurable peaks and 
three null results, all of which can be used to constrain 
the magnetic structure. The orbital contributions to the 
magnetism in CeFe 2 are negligible compared to the spin 
contribution for both Ce and Fe |S]. The non-resonant 
magnetic diffraction cross section is thus dominated by 
the projection of the spin onto the direction (S2) perpen- 
dicular to the scattering plane. The measured quantity 
is the ratio of diffraction intensity from the antiferromag- 
netic and the lattice Bragg peaks, expressed as 



j I h k I \ 

lLAr(h,k,l) 



m e c 2 8 / e e iqr 



sin 2 (26>). (1) 



Here (h, k, 1) are Miller indices of the unit cell before 
antifcrromagnetic period doubling. The summations in 
the numerator and denominator run over all scattering 
sites in the magnetic and lattice unit cells, respectively, 
and the factor of eight accounts for the difference in size 
between the magnetic and lattice unit cells. The x-ray 
energy is hu, while 20 is the diffraction angle of (h/2, 
k/2, 1/2), and f m and f e are the magnetic and atomic 
form factors |S]. 

In order to constrain the magnetic structure we re- 
sort to a general treatment of antifcrromagnetic order on 
a face centered lattice [22]. The ordering doubles the 
unit cell along all three axes in real space. For a face 
centered lattice the basis therefore increases from one to 
eight points, forming a bipartite lattice consisting of basis 
points 1-4 and l'-4', with the condition that the primed 



and unprimed points be magnetically distinguishable. In 
Fig. [4] we choose the assignment of the basis points that 
is most natural to a structure layered along (111). Bear 
in mind that the Laves basis consists of six atomic sites 
(four Fe and two Ce) at each point of the bipartite lattice. 

The antiferromagnetic structure proposed for 
Ce(Fei_ a; Co a: )2 on the basis of neutron and reso- 
nant x-ray diffraction measurements [51 [TJ, consists 
of ferromagnetic kagome sheets of 3e-Fe spins (each 
carrying 1.61 /iB) polarized along (111), with the spins 
inverted on adjacent sheets. The Ce spins (0.13 /iB) 
are also parallel to (111) and are ordered antiferro- 
magnetically within each (111) plane. The lb-Fe spins 
(1.12 /iB) are polarized within the (111) plane. This 
can be understood as the result of frustration, arising 
from the ferromagnetic coupling of the lb-Fe spins to 
the two oppositely polarized neighboring kagome sheets. 
Considering only the coupling of the lb-Fe spins with 
their (Ce and 3e-Fe) nearest neighbors, all orientations of 
the spins within the plane are degenerate. In the model 
of Ref. [7J the sheets of lb-Fe spins are fcrromagnetically 
polarized in the plane, so that the lattice sites 1-4 have 
identical spin orientations, which are opposite to the 
spins on sites l'-4'. However, it is straightforward to 
show that any such model gives a magnetic structure 
factor of zero for all six of the observed reflections 
plotted in Fig. [3j 

The model proposed for Ce(Fe 1 _ a; Co 2 ;)2 is therefore 
not directly transferrable to antiferromagnetic CeFe2 un- 
der pressure. However, this model can conform to our 
data with one revision. The fact that the Fe spins are 
not all co- linear is established [7J, and the assumption 
that the lb-Fe spin is misaligned with the 3e-Fe spins is 
consistent with the rhombohedral distortion. However, 
the interaction between lb-Fe spins and their orienta- 
tion within the (111) plane have not been constrained. 
If this interaction is antiferromagnetic rather than ferro- 
magnetic, then the triangular plaquettes of lb-Fe sites 
will be magnetically frustrated, the bipartite lattice sites 
1-4 will not have identical spin orientations, and the ob- 
served reflections acquire non-zero structure factors. 

As explained in detail in the Supplemental Material, 
this leads to a model (Fig. |4| that is consistent with 
our data. The planes of lb-Fe sites consist of antifer- 
romagnetically aligned pairs and two spin polarization 
axes separated by £ = 60°. There is degeneracy associ- 
ated with the choice of pairs, giving rise to a magnetic 
domain structure. Each domain contributes to some but 
not all of the reflections shown in Fig. [3]d and for simplic- 
ity we assume that the three domain types are equally 
populated. We further assume that the angle cp in Fig. [1J; 
differs by 120° between domains. The diffraction intensi- 
ties calculated from this model are in decent quantitative 
agreement with our data (Fig. ^p), including the three 
null results, and allow us to determine the lb-Fe mo- 
ment to be 1.80±0.04/d3. The agreement between model 
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FIG. 4: (color online). Antiferromagnetic structure of CeFe2- 
(a) Schematic showing the bipartite sub-lattices (1-4) and (1'- 
4') for a face centered antiferromagnet [22]. Each sub-lattice 
point is associated with a complete Laves basis, which is ex- 
plicitly drawn on sub-lattice point 1. The spin orientations 
on site 1' are the inverse of those on site 1, and likewise for 
the other pairs, (b) The 3e-Fe spins on each kagome sheet are 
ferromagnetically aligned, (c) The lb-Fe spins form a plane of 
frustrated triangular plaquettes. The degeneracy associated 
with the choice of antiferromagnetic pairs leads to a magnetic 
domain structure; only one domain is drawn here as an ex- 
ample, (d) Monte Carlo simulation of the generic magnetic 
phase diagram, which is controlled by competition between 
the antiferromagnetic Ce-Ce exchange interaction (J3), and 
the ferromagnetic Fe-Fe interaction (J2). The spin structure 
for CeFe2 in both the ferromagnetic and antiferromagnetic 
phases is drawn for the lb-Fe site and its twelve nearest neigh- 
bors. 



and measurements can be further improved by assuming 
unequal magnetic domain populations and/or by opti- 
mizing the angles ip and £, but we feel such an analysis 
is unjustified within our limited data set. 

We can now understand the phase diagram by consid- 
ering the magnetic interactions. In both the high and 
low pressure phases, the Ce spins are anti-aligned with 
the majority of the Fe spins 6 8, 23J . This suggests that 
an antiferromagnetic coupling (with exchange energy J\ ) 
between the Ce and Fe spins is the dominant magnetic 
interaction. This coupling may be expected to arise from 
the double exchange interaction between the localized 
Fe moments and the itinerant Ce electrons [23], and is 



known to exist in Laves systems [5[ [35] . The bonds be- 
tween Fe nearest neighbors are ferromagnetic (exchange 
energy J2). Finally, the Ce spins are antiferromagneti- 
cally aligned in the high pressure phase. The proximity of 
the magnetic phase transition in CeFe2 to ambient pres- 
sure, and the insensitivity to the interatomic distances 
(Fig. [I]) , are strong indications that the magnetic inter- 
actions will not be qualitatively different on opposite sides 
of the transition. We therefore propose that the Ce-Ce 
interaction is always antiferromagnetic (exchange energy 
J3), and that the shifting balance between the antifer- 
romagnetism of the Ce spins and the ferromagnetism of 
the Fe spins drives the transition. 

To test this proposal we construct a semiclassical 
model, taking into account only the nearest neighbor in- 
teractions Ji, J 2 and J3. We numerically determine the 
phase diagram both within a mean field approximation 
and using a classical Monte Carlo simulation (see Supple- 
mentary Material) , assuming that the Ce and 3e-Fe spins 
will be aligned along (111). The resulting phase diagram 
is shown in Fig. [4}i. Even at this level of approximation, 
the theoretical model clearly reproduces the qualitative 
features of the experimental phase diagram. 

The generic magnetic phase diagram in Fig. [T] is now 
explained. The low temperature, ambient pressure phase 
of CeFe2 is ferromagnetic, but it sits close to an antifer- 
romagnetic instability and J3/J2 is close to the critical 
value. Since Al, Co, Ru and Ir all replace Fe upon be- 
ing doped into CeFe 2 , they affect the Fe-Fe bonds more 
strongly than the Ce-Ce bonds. The dopants which end 
up on Fe-lb sites decrease the energetic cost of the an- 
tiferromagnetic stacking of Fe-3e layers and thus effec- 
tively increase J3/ ' Ji- On the other hand, the application 
of pressure mainly affects the interatomic distances, and 
therefore affects both the ferromagnetic and the antifer- 
romagnetic bonds. The latter is more sensitive to the 
changing overlap integrals due to the itinerant character 
of the Ce electrons, thus shifting J3/J2 towards higher- 
values under pressure. The horizontal scaling of various 
magnetic phase diagrams for chemical doping and ap- 
plied pressure results from the fact that there is only a 
single parameter J3/J2 which controls the transition be- 
tween the ferromagnetic and antiferromagnetic phases. 
By contrast, the thermal transition into the paramag- 
netic phase is controlled by J\. This strong Ce-Fe spin 
interaction is little affected by either pressure or doping, 
and the vertical axes of the different phase diagrams thus 
line up without any additional scaling. 

The picture which emerges for CeFe2 contains multiple 
magnetic energy scales which compete in a landscape of 
geometrical frustration, and the magnetic structure is de- 
termined by the balance between the weakest two interac- 
tions. This solution holds for both chemical substitution 
and applied pressure and may serve more generally for 
other geometrically frustrated pyrochlore magnets. The 
phase transition occurs when the effective interaction be- 
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tween kagome layers turns from ferromagnetic and medi- 
ated by the Fe spins, to antiferromagnetic and mediated 
by the Ce spins. In the ferromagnetic phase, the Fe-Fe 
interactions are satisfied and the Ce-Ce interactions are 
frustrated. The rhombohcdral distortion tilts the play- 
ing field in favor of the Ce-Ce interactions by isolating 
the lb-Fe spins, with the result that in the antiferromag- 
netic phase the Fe-Fe interactions are partially frustrated 
while the Ce-Ce interactions are mostly satisfied. Under- 
pinning the whole phase diagram is the dominant antifer- 
romagnetic Ce-Fe interaction, which depends on strong 
hybridization between the itinerant Ce 5d and the more 
localized Fe 3d orbitals. 
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METHODS 

Single crystals of CeFe2 were grown from a Ce-rich bi- 
nary melt with initial composition Ceo.gFeo.4. High pu- 
rity Ce (Ames Lab) and Fe were sealed into a three-cap 
Ta crucible [10] and subsequently sealed into a silica am- 
pule. The ampule was heated from room temperature to 
1100 °C over 6 hours, cooled to 950 °C over 3 hours and 
then slowly cooled to 700 °C over 120 hours. Once at 700 
°C the excess liquid was decanted from the single crys- 
tals |10| . X-ray diffraction measurements were performed 
at beamline 4-ID-D of the Advanced Photon Source. A 
double-bounce Si (111) monochromator and a pair of Pd 
mirrors produced a focused beam of 20 keV x-rays free 
from contamination by higher harmonics. CeFe2 crystals 
of typical dimensions 70 x 70 x 40^m 3 were loaded in 
the diamond anvil cell (DAC) in an argon atmosphere 
to avoid oxidation. The pressure medium was a 4:1 
(volume) methanokethanol mixture. A polycrystalline 
silver grain was used as an in-situ manometer |19j . A 
helium-mcmbranc-tuncd DAC was used to continuously 
vary pressure at the cryostat base temperature of 3.5 K. 
X-ray scattering was carried out in the transmission ge- 
ometry and within the vertical scattering plane for high 
momentum space resolution. The use of 70° full cone 
Boehler-Almax diamond anvils allowed access to a wide 
range of reciprocal space. A total of four samples were 



studied under pressure at 3.5 K. The FWHM of the sam- 
ple rocking curve never exceeded 0.10° under pressure. 
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